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ABSTRACT 

Lyman Alpha Emitters (LAEs) are usually thought to be powered by star formation. 
It has been recently reported that a fraction of LAEs at redshift z ^ 3 — 4 hosts an 
Active Galactic Nuclei (AGN). If an AGN is present it could be obscured and unde- 
tectable in X-rays, but yet dominate the Lya luminosity. We examine the properties 
of these AGN-powered LAEs at high redshift {z ^ 6) using radiative transfer cosmo- 
logical simulations and obtain a reliable criterion to identify them from their observed 
Lya line and surface brightness. We find that these sources should have: (a) negative 
line weighted skewness, < 0, and (b) surface brightness profiles FWHM ^ 1.5". 
This parameter space cannot be populated by starburst LAEs. Thus, LAEs satisfying 
this criterion would be strong candidates for the presence of a hidden AGN powering 
their luminosity. 

Key words: galaxies :high-redshift - fundamental parameters - evolution - abundances 
- stellar content 



1 INTRODUCTION 

Lyman Alpha Emitters (LAEs) are galaxies emitting 
a prominent Lya 12 16 A emission line. Narrow band 
imaging techniques have allowed the discovery of thou- 
sands of LAE candidates in an extended redshift ran^e 



2-7 ( Cowie & Hu, Il998l : 



20041: IShimasaku et al.1 l200, 



Malhotra k Rtioadd 
Ouchi et al.1 l2005l . l2008l : 



Kashikawa et al.l 1 20061 ) and hundreds of them have been 



spectroscopically confirmed (IHu et al.l [20o3; IPawson et al.l 
l2004l : lOuchi et al.ll2008l : iKashikawa e"taDl2006l ). 

While the most natural explanation for the large Lya 
luminosities is associated with a star forming activity, inter- 
estingly, a fraction of LAEs seem to host an Active Galac- 
tic Nucleus (AGN). This important property is now con - 
firmed by X-ray detectio ns report e d by Gawiser et al.l (|2006l l 
at z - 3, and by .Zheng et al.1 (l2010l) at z ^ 4 5 . The 
source J033127.2-274247 discovered by IZheng et all (|2010l ) 
shows a strong soft band (0.5 — 2 keV) X-ray flux with 
Lx = 4.2 X 10^^ ergs~^, and has been spectros c opica lly con- 
firmed as an unobscured AGN. lOuchi et al.1 (|2008l ) found 
broad emission lines - a typical AGN signature - among 
photometrically selected LAEs at z = 3.1 and z = 3.7. 
The fraction of LAEs showing broad emission line is small 

1%); however, the brightest LAEs (L^ ^ lO^^ '^erg s"^) 
always appear to host an AGN. 

Unfortunately, at high redshifts, where assessing the 
presence of a central accreting black hole would be of out- 
most importance, e.g. to study the precursors of the super- 
massive black holes. X-ray fiuxes are too faint to be used as 
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a discriminating tool. In addition, it is likely that in the ma- 
jority of these sources X-ray emission from the central black 
hole is obscured by a dense absorbing gas layer. One would 
then like to be able to infer the presence of an AGN directly 
from the properties of the Lya emission, which most often is 
the only information available for these distant sources. T he 
giant LAE at z = 6.6 named Himiko (jOuchi et al.ll2009l ) is 
suspected to contain a hidden AGN due to its prominent 
size and Lya luminosity. The lack of detection in either X- 
ray, MIR, sub-mm and radio bands are not sufiicient to dis- 
card such possibility. In this work, based on high resolution, 
radiative transfer (RT) cosmological simulations of a LAE 
similar to Himiko, we derive a criterion to unambiguously 
assess if a generic LAE is powered by an AGN or by a star- 
burst. 



2 NUMERICAL SIMULATIONS 

We start by running a cosmol ogical SPH hyd rodynamic sim- 
ulatiorQ using GADGET-2 (ISpringell l2005l ) and extract a 
snapshot at z = 6.6. We simulate a (10/i~^Mpc)^ volume 
with 2 X 512"^ baryonic+dark matter particles, giving a mass 
resolution of (1.32,6.68) x 10^ M© for (baryons, dark mat- 
ter). We post-process the snapshot by runn i ng UV and X- 
ray RT using LICORICE jBaek et al.ll2009l . boiol ). so that 
the remaining average H i fraction, xhi, in the intergalac- 
tic medium (IGM) is 0.1, in agreement with current data 



We use the recent WMAP7+BAO+ifo cosmological param- 
eters: Qm = 0.272, Qa = 0.7 28, Qt = 0.0455, h = 0.704, 
erg 



: 0.807 dKomatsu et al.ll201ll V 
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Figure 1. Gas density field in a slice of thickness 1 kpc and 
proper size 120 kpc, centered on the most massive halo with = 
1.17 X 10^^ Mq at z = 6.6. The colorbar shows the density values 
in [cm""^]. 

(jMortlock et al.l [20111 ). By using a friends-of- friends algo- 
rithm we identify the most massive halo, of total mass 
Mh = 1.17 X 10^^ Mo, gas mass Mg = 2 x lO^^M©, and 
virial radius rvir ~ 20 proper kpc (pkpc). Fig [1] shows the 
gas density map of the halo, connected to the filamentary 
dense structure. Then we select a (0.625 h~^Mpc)^ comoving 
volume centered on the halo and re-simulate the UV/X-ray 
RT at higher spatial resol ution on an adap tive grid with 
Nmax = 8 according to the iBaek et al] (| 20091 ) schem^l, giv- 
ing a minimum RT cell size of 0.114 pkpc. 

We have run simulations for three different Spectral En- 
ergy Distributions (SED) of the central source, namely (a) 
starburst, (b) Compton-thick, (c) Compton thin AGN types, 
keeping all other conditions the same, and setting the ioniz- 
ing photons escape fraction, fesc = 0. For the starburst case, 
radiation is emitted isotropically with a bolometric luminos - 
ity obtained from STARBURST99 (|Leitherer et al.1 Il999l ) 
assuming a continuous star formation rate of lOMoyr"^, 
metallicity Z = 10~^, a Salpeter initial mass function a = 
2.35 in the mass range 1 — 100 M©. The resulting photoion- 
ization rate is Q = 3 x 10^^ s~^. For AGNs we use t he com- 
posite spectra from radio to X-rays following Shang et al.l 
(|201lh , with an X-ray upper cut-ofl0 at 0.5 keV; the spec- 
trum is normalized so that all cases have the above pho- 
toionization rate Q. Compton-thick AGNs emit their ioniz- 
ing radiation (hu > 13.6 eV) in a cone with opening angle 
= 45° with respect to the polar axis perpendicular to the 
line of sight; for Compton-thin AGNs, in addition, energetic 
X-ray photons {hu > 0.1 keV) are allowed to emerge from 

^ The Nmax = 8 value ensures that the RT cell contains only few 
particles: of the total number of cells, a fraction (0.31, 0.28,0.16) 
contain (1,2,3) particles. The number of UV photons in a photon 
packet has to be smaller than the number of H atoms in the cell. 
Even with adaptive integration time steps, this is computationally 
too demanding. Therefore gas particles outside the re-simulated 
box are set to have xiii=0.1 and T = 10^ K for simplicity. 
^ More energetic photons have a mean free path larger than the 
galaxy size and hence do not contribute to Lya emissivity 



the dusty torus (size ^ 1 pc) and propagate isotropically 
into the ISM. The rest-frame AGN X-ray luminositjO in the 
energy band corresponding to the 0.5-2 keV observed band 
is 2.42 X 10^^ ergs"^ 

In addition to H recombinations, Lya photons can be 
produced also by X-rays via de-excitations. Soft X-rays (0.1- 
0.5 keV) have a relatively short (1 — 100 kpc) mean free path 
in the interstellar medium of z ^ 6 galaxies. Once the X-ray 
photon is absorbed by an H atom, the primary electron of 
energy Ep = hiy — 13.6 eV can produce further (secondary) 
ionizations, excitations and heating of the gas; de-excitations 
from upper energy levels to the ground state eventually re- 
lease Lya photons. The fractiona l energy of the prim ary elec- 
tron that goes into Lya photons (jValdes FerraralfioOS ) de- 
pends on xhi: if xm> 0.99, f^{xm) ~ 30%; such efficiency 
decreases to (xhi) = 0.01 for xm= 0.1. 

Following X-ray photon absorption we compute the ad- 
ditional production of Lya photons by de-excitations as 

AiVLyc. = Nabs^faixmi), (1) 

where Nabs is the absorbed number of photons and hua 
the Lya photon energy. During UV/X-ray RT, we store 
the number of recombination and de-excitation Lya pho- 
tons produced by each SPH particle. 

Finally Lya RT is ru n using the LICORICE Lya mod- 
ule (ISemehn et al.1 l200/t ) on a volume of 0.625 x 0.625 x 
10/i-^Mpc^ The size in the sky plane is sufficiently large 
(20" X 20") to model the surface brightness (SB) map of 
a LAE; the larger size (240 bins with rest frame size 0.1 
A, centered on Aa) along the line of sight /frequency axis 
is required to properly account for the Lya IGM transmis- 
sivity. We interpolate all gas physical properties on a fixed 
128 X 128 X 2048 grid using the SPH kernel, thus obtaining 
a spatial resolution of 0.918 pkpc. This is almost identi- 
cal to the p ixel scale achieved by the Subaru/Suprime-Cam 
(iMivazaki et al. 2002) in the Subaru/XMM-Newton Deep 
Sky Survey |Ouchi et al.ll2008l V We compute the Lya lumi- 
nosity from each cell, L^^\ by summing all Lya photons 
emitted by the SPH particles in the cell, 

= ^ E ^Ly., (2) 

where Nceii is the number of SPH particles in the cell and i 
is the particle index. If L^^^ ^ lO^^ergs"^, the cell is con- 
sidered as a Lya sourc^l, and we cast Lya photon packets 
from it. 



3 RESULTS 

We show the total (recombination + deexcitation) Lya lumi- 
nosity. La, of the simulation box as a function of time in Fig. 
[2] For the starburst case La increases rapidly as the gas gets 
progressively ionized by the central source, approaching an 
asymptotic value of 5 x 10^^ erg s~^, set by the equilibrium 
between recombination and photoionization rates. Breakout 

^ These sources are then 10 times fainter than the detection limit 
of XMM-Newton X-ray in the same band. 

^ Even for a detection limit as low as 10~^-'^erg s~-'^cm~^arcsec~^, 
cells with L*^^^ < lO^^ergs"-*^ contribute negligibly. 



LAEs powered by AGN 3 




0.4 0.6 0.8 
Time [Myr] 



1.2 



Figure 2. Time evolution of the total Lya luminosity of the sim- 
ulation box. The various curves refer to starburst (red dashed 
line), Compton-thin (black solid) and Compton-thick (black dot- 
ted) SED cases. 



of the ionization front (IF, defined by the position where as 
xmi = 0.5) from the galaxy occurs at 0.25 Myr. 

For the Compton-thick AGN the evolution is different. 
Initially UV photons are absorbed locally, but X-rays pene- 
trate to larger distances thus making the IF much smoother 
than for the starburst SED; X-rays absorbed beyond the 
IF effectively produce additional Lya photons. The IF also 
propagates faster than in the starburst case, as the same 
photoionizing rate Q is now confined in a cone. While the 
IF travels in the ISM, X-rays boost the Lya luminosity by 
several times with respect to the starburst SED; however, 
beyond the halo virial radius the low IGM density strongly 
reduces the Lya emission; hence after 0.1 Myr the total Lya 
luminosity drops as the ionized volume where recombina- 
tions occur is confined to a cone. Thus we can identify two 
separate regimes for Lya luminosity, namely (a) an X-ray 
dominated regime (IF in the ISM), and (b) a recombination 
dominated regime (IF in the IGM). 

For the Compton-thin AGN, in which X-rays are emit- 
ted isotropically, de-excitations in the neutral ISM outside 
the torus can largely amplify the Lya luminosity with re- 
spect to the Compton-thick case as a larger fractional energy 
of X-rays can be converted into Lya photons. Therefore, the 
case of Compton-thin AGN produces the highest luminosity 
among all cases, and converges slowly than the Compton- 
thick AGN case. 

X-rays produce Lya very efficiently in the initial phases 
of photoionization, when the ISM is largely neutral. This can 
be appreciated from Fig.O showing ionization fraction, Lya 
luminosity, and SB maps extracted 0.06 Myr after source 
turn-on for the three source cases. The snapshot at 0.06 
Myr shows the most distinctive feature between starburst 
powered and AGN powered LAEs. The starburst SED pro- 
duces a very sharp IF as expected, while the AGN SEDs 
are characterized by smoother IFs. For the Compton-thin 
case, the two outer contours (xhi= 0.99, 0.999) are quite 
isotropic; nevertheless, the cone-shaped HII region is clearly 
visible. 



In the starburst case, Lya photons are produced only 
inside the HII region through recombinations. As the re- 
combination rate is proportional to n^, the Lya luminos- 
ity map closely resembles the density distribution, i.e. it is 
brighter at the center; this property also holds for the two 
AGN cases within the cone. Apart from these dense, ionized 
regions where recombination-driven Lya production domi- 
nates, Lya are created by X-rays also near and beyond the 
IF. As Compton-thin AGNs emit X-rays (hu > 100 eV) 
isotropically, a copious amount of Lya photons comes from 
regions shielded from the UV flux. 

As for the SB distribution (Fig.[3l bottom row), AGNs 
show longitudinally extended shapes, while the starburst 
case is more isotropic. Even though Lya emission is stronger 
in the central regions of the galaxy, SB maps do not ex- 
actly match the corresponding La maps. This is due to 
the fact that most Lya photons experience considerable fre- 
quency and spatial shifts through repeated scatterings with 
H I atoms in the ISM. After being scattered ~ 10^ times, the 
photon frequency is shifted away from the Lya line center, 
allowing escape even if the H i column density is as high 
as A^Hi = lO^^cm"^. None of pixels in the starburst and 
Compton-thick AGN cases is above the current detection 
limit of narrow band images, 10~^^ ergs~^ cm~^ arcsec"^; 
for the Compton-thin AGN instead, the region within the 
central 1.5 arcsec^ can be detected. The halo is connected 
with dense filamentary structure as shown in Fig[T] but we 
checked that the total amount of HI in the filament is not 
enough to produce detectable Lya emission. 



3.1 Skewness evolution 

From the observed emission Lya line profile we can compute 
the line skewness, 



■xffi.. 



(3) 



where fi is a 2D of the flux, Xi is the pixel coordinate, 
/ = fi^ cind X, a are the mean and dispersion of Xi. 
Due to its higher sensitivity to profile asymmetries, the 
weighted skewness, Sw = S{Xio,r — Aio,b), where Aio,r and 
Aio,b are the wavelengths where the flux drops to 10% of the 
emission peak va l ue on the red and blue sides of the line 
(|Kashikawa et all l2006l ) . is more often used. Typically the 
Lya line from high- 2; sources, differently from other emis- 
sion lines such as [O II] or [O III] at lower redshifts, has 
Sw > 0. Thus, Sw can be an useful tracer of high- 2; objects. 

For each SED case, we compute Sw for 5 snapshots 
(t = 0.02,0.04,0.06,0.08,1 Myr), the first four (last) cor- 
responding to the X-ray (recombination) dominated regime. 
It is well known that for a spherical, static case (iNeufeldl 
I1990I : iDiikstra et al1l2006l ) , the emerging line profile shows a 
double peak shap43, with the peaks getting further apart as 
A^Hi is increased. The blue wing is then suppressed as pho- 
tons propagate through the IGM if the source is located at 
sufficiently high redshift z ^ 5; the red wing can be instead 
transmitted to the observer, giving Sw > 0. 



^ The size of the core is x = {h> — h'oi)/^D ~ 3 where ujj = 
(vth/c)iya is the Doppler frequency shift and Vth is thermal ve- 



locity dispersion of the gas. 



4 S. Baek, A. Ferrara 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 



# 





37.4 37.8 38.2 38.6 39 39.4 39.8 40.2 40.6 37.4 37.8 38.2 38.6 39 39.4 39.8 40.2 40.6 37.4 37.8 38.2 38.6 39 39.4 39.8 40.2 40. 




Figure 3. Maps of the H l fraction with contour levels xhi^O.Q, 0.99, 0.999 (top row), Lya emissivity (middle), and SB with contour 
levels 10-21+^', j = 0, ..,3ergs-i cm arcsec 2, smoothed with a gaussian of FWHM=0.8" (bottom) for the (a) starburst (left column), 
(b) Compton-thick AGN (middle), and (c) Compton-thin AGN (right) cases. Snapshots are taken 0.06 Myr after source turn-on. Maps 
have size 60x60 pkpc and are 3 pkpc-thick through the center of the most massive halo in the box. 



During the X-ray dominated stage of Lya production 
A^Hi is very high (10^^~^°cm~^), corresponding to an optical 
depth at line center = 10^~^. The peak of the red wing is 
shifted farther to the red and the profile decline on the blue 
part of the peak becomes less steep than on the red part, 
thus decreasing Sw • Although our simulated gas density and 
velocity field are far more complex that those underlying the 
Neufeld solution, the above interpretation continues to hold. 
For this reason we find small or even negative values of Sw 
As photoionization proceeds, A^hi drops below lO^^cm"^, 
the peak position gets closer to the line center, and the pro- 
file becomes highly skewed as typically observed in LAEs. 
Most cases have Sw < initially, later increasing as A^hi 
decreases. 



3.2 Surface brightness evolution 

Star- forming galaxies are more compact at high z: typ- 
i cally, at z ^ 6 they have a half-light radius of 1 kpc 
teouwens et al Therefore the Lya emitting area is 

smaller than the angular resolution ( ^ 0.8") of narrow band 



images and the SB radial profile from the brightest pixel ap- 
proximates a Gaussian point spread function (PSF) of s.d. 
a = seeing/2.35. If Lya emission is powered by X-rays, how- 
ever, it can be spatially more extended, as already stated. 



Fig.|4]shows the simulated SB profiles for the three SED 
cases. For the starburst, the SB profile is initially relatively 
broad (FWHM ^ 1.2") since Lya photons are very scattered 
in space; later on, the FWHM decreases and approaches the 
resolution scale expected for a Lya point source; even in 
the early evolutionary phases, though, the FWHM never 
exceeds 1.4". On the contrary, for both AGN cases several 
pixels are brighter than the detection limit when the Lya 
photon production is dominated by X-rays (Fig. |4|. The 
width is ^ 1.5" for the Compton-thin AGN case. After 1 
Myr, though, the flux fades below the detection limit and 
the profile becomes narrower since the total Lya luminosity 
decreases. 
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Figure 4. SB profile evolution through the brightest pixel of the maps in Fig.[3]at t = 0.02, 0.04, 0.06, 0.08, 1.00 Myr 
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Figure 5. FWHM of SB profile vs. skewness. From left to right, 
each model has 5 data points at t = 0.02, 0.04, 0.06, 0.08, 1.00 
Myr. The shaded area denotes AGN-powered LAE identification 
criterion Sw < and FWHM ^1.5", 



3.3 Identification criterion 

From the previous results we can build a robust criterion 
to identify AGN-powered LAEs, as shown in Fig. [5] where 
we plot the evolutionary tracks of the three SEDs in the 
FWHM vs. Sw plane. As Sw is anti-correlated with A^hi, it 
increases with time for all models. Thus Sw < alone cannot 
be used to uniquely identify AGN-powered LAEs. However, 
the starburst FWHM cannot exceed 1.4", independently of 
the assumed SFR as long as it is sufficient to allow detection 
of the source as a LAE (La > lO"^^'^ ergs~^). Therefore, the 
criterion defined hy Sw < and FWHM ^ 1.5", can be 
safely used to identify the class of LAEs powered by AGN, 
if they exist. 



4 CONCLUSIONS 

The presence of nuclear black holes, progenitors of the super- 
massive ones discovered at 2; = 6 — 7, in high redshift galaxies 
is expected on several grounds. Yet direct detection of their 
X-ray emission has so far proven very challenging, partic- 
ularly at 2; > 6. In addition some of these AGNs can be 
heavily obscured and show no broad emission lines. Thus it 



would be useful to be able to identify AGN-powered LAEs 
purely from their observed Lya properties. 

We have explored this possibility using RT cosmolog- 
ical simulations of one of the most luminous LAE known, 
Himiko. From the simulated Lya luminosities, line profile 
and surface brightness we have been able to isolate key dif- 
ferences among three plausible SEDs of the central source 
(starburst, Compton-thick or Compton-thin AGN types) 
with a fixed photoionization rate Q = 3 x 10^^ s~^. 

From the results we have built a robust criterion to iden- 
tify AGN-powered LAEs. We find that these sources should 
have (a): negative line weighted skewness, Sw < 0, and (b) 
surface brightness profiles FWHM ^ 1.5"; this parameter 
space cannot be populated by starburst LAEs. Thus, LAEs 
satisfying this criterion would be strong candidates for the 
presence of a hidden AGN powering their luminosity. Note 
that this criterion purely depends on the observed properties 
of the Lya line and does not require additional information 
from, say, UV continuum or X-ray data. According to such 
criterion, Himiko, due to its high skewness, Sw = 13.2 and 
somewhat too narrow SB profile (Fig|5]), is not predicted to 
host an AGN. 

Clearly the success of the proposed strategy depends 
on the sensitivity of the observations. A sufficiently deep SB 
map will allow to accurately measure the SB radial profile 
over a wide range of distances where the differences between 
starburst and AGNs become more appreciable (see Fig. 2]). 
In this sense, the advent of the new generation of instru- 
ments, as e.g. MUSE, will be crucial. 
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